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Nonextensive quantum gas distributions are investigated on the basis of the factorization hy- 
pothesis of compound probability required by thermodynamic equilibrium. It is shown that the 
formalisms of Tsallis nonextensive statistical mechanics with normalized average give distribution 
functions for standard bosons and fermions obeying Pauli principle. The formalism with unnormal- 
ized average leads to a intermediate quantum distribution comparable to that of fractional exclusion 
statistics, with Fermi surface at T = depending on the parameter q. 
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I. INTRODUCTION 



It is believed that long-range correlations can often (if not always) be related to fractal or chaotic space-time 
structures and also to power laws and nonextensive phenomena Jl]-^]. That is why the fractal inspired nonextensive 
statistical mechanics (NSM) Q is receiving much attention and widely considered as a possible valid theory to 
describe systems with complex correlations. Considering the fractal or chaotic behaviors and nonextensive or fractional 
effects in quantum systems like, for example, correlated electrons in superconductor, solitons and quasi-particles in 
condensed matters we naturally expect a quantum statistical mechanics within NSM. One of the first steps 

in this direction is taken by Buyiikkilig et al jl2| in giving the generalized quantum distributions with the so called 
factorization approximation. 

The derivation of the conventional quantum one-body distribution from the compound many-body distribution 
is straightforward within Boltzmann-Gibbs statistics because the latter considers only short range interactions with 
additive entropy and energy. This allows the many-body density operator p to be factorized into the product of 
one-body density operators pi of N particles of the system : 

N 

P = \{pi- (1) 

i=l 

But in nonextensive statistical mechanics, this kind of factorization of compound probability is impossible if we still 
consider, for independent particles or ideal gas, extensive energy with total hamiltonian H given by 

N N 2 

= (2) 

1=1 1=1 

2 

where Hi — + Vi is the one-body hamiltonian of i th particle, pi its momentum and Vi its potential energy. 
To overcome this difficulty, a so called factorization approximation is proposed which imposes Eq.(|l|) in keeping 
Eq. (||) ■ This approximation allows to obtain a nonextensive quantum distribution of NSM with complete distribution 
(i.e. Trp = 1) and unnormalized average (i.e. x = Trp q x where x is certain operator Q). It is, explicitly or not, 
employed in most of the applications of NSM to the cases like, among others, boson and fermion systems [p^|-^4|, 
the polytropic model of galaxies, solar neutrinos, peculiar velocity of galaxy clusters, electron plasma turbulence 
(for updated comments on these works, we refer to reference [|l5| ), and the application of nonextensive blackbody 
distribution to laser physics |lqj . Although these applications clearly evidence the existence and the important role 
of the ^-exponential distribution in the nature, the justification of the passage from many body g-distribution to 
one-body one is missing. Recently, a direct analysis |l7j shows that, assuming Eq.(|l|) and Eq.(||) at the same time 
means neglecting a correlation energy which is in general not negligible. A detailed computation for A-oscillator 
system [ [l8[ also shows that the partition function given on the basis of Eq. ([!]) is completely different from that given 
by using Eq.(|). 

In this paper, the factorization hypothesis is discussed from a new point of view relative to Abe's general pseu- 
doadditivity [ pl| for nonextensive systems in thermal equilibrium. This discussion leads to the conclusion that Eq.(Q) 
must be respected in all coherent and exact treatments of equilibrium systems. On this basis, quantum distributions 
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are derived within different formalisms of NSM : 1) complete distribution (Trp = 1) with the normalized average 
x = Jrpx [Q; 2) incomplete distribution with g-normalization Trp 9 = 1 and the normalized average x = Trp q x [fj"7| ; 
3) complete distribution with unnormalized average x = Trp q x 0. 



II. THERMAL EQUILIBRIUM AND FACTORIZATION OF COMPOUND PROBABILITY 

In Boltzmann-Gibbs statistics, Eq.(|l|) is a natural result of the short-range interactions or of the independence of 
the subsystems which necessarily lead to extensive entropy S and internal energy E. With the temperature definition 
l/T = (3 = J|t (Boltzmann constant k=l), the extensivity of S and E allows to verify the zeroth law or the existence 
of thermal equilibrium. So, the uniqueness of Gibbs-Shannon entropy S = — YliPi m Pi subject to Shannon-Khinchin 
axioms jl9| can be regarded from a different viewpoint. We can also state that, given Gibbs-Shannon entropy, 
the existence of thermal equilibrium necessarily leads to the additivity of entropy or to Eq.(|]J). This statement 
seems self-evident and useless for Boltzmann-Gibbs statistical mechanics, because Eq.(Q) and the axiom about the 
additivity of entropy are obvious in view of the independence of subsystems. But it is not the case in NSM . The 
pseudoadditivity axiom S — S\ + S2 + (1 — q)SiS2 @ for a composite system containing two subsystems 1 and 2 
with respectively entropies 6*1 and S2 , is in fact a consequence of Eq. (|l|) supposed Tsallis entropy S = — 1 ~ T y" . As a 
consequence, the nonextensive theory is entirely subject to the "independence" of subsystems. The problem is that, 
in nonextensive physics, the subsystems are in general dependent on each other. So Eq. ([!]) or the pseudoadditivity is 
in this case somewhat artificial and seems to hide something which may restrict the validity of the theory. In many 
applications of NSM, Eq.fll]) is sometimes respected, sometimes not, sometimes combined with the incompatible 
extensive hamiltonian Eq.(Q) |l2|,|l8|,^0| . This confused situation is in our opinion due to the fact that Eq.(jl|) is not 
clearly founded for nonextensive statistics subject to Tsallis entropy. 

Very recently, Abe |2lf| proposed a general pseudoadditivity for entropy required by the existence of thermal 
equilibrium in composite nonextensive systems : f(S) = f(Si) + f(S2) + Xf(Si)f(S2) where / is certain differentiable 
function satisfying /(0) =0 and A a constant depending on the nature of the system of interest. So for a system 
containing N subsystems, the thermal equilibrium requires following additivity : 

N 

ln[l + A/(S)] = £>[l + A/(#)]- (3) 

i=l 

On the other hand, Eq.(||) applied to Tsallis entropy means f(S) = S and A = 1 — q plj , which directly leads 
to lnTrp 9 = ^ i=1 hiTrp| ° r Eq. ([!]). This interesting result raises Eq.(Q) from the level of a special assumption 
for "statistically independent" subsystems to the level a general theorem of equilibrium thermodynamics with cither 
Boltzmann-Gibbs or Tsallis entropy. We does not need any more the independence of subsystems to write Eq. (Q) which 
must be respected by all probabilities of the systems (extensive or nonextensive) that can have thermal equilibrium. 
Consequently, all calculations based on the factorization theorem or using one-body g-exponential distribution are, 
as a matter of fact, correct and exact applications of NSM. And all calculations of NSM based on Eq.(|J) [not 
compatible with Eq.(Q)] should now be regarded as "extensive energy approximation" which may be very different 
from the exact treatment and should be employed with great care. 

In what follows, we discuss the quantum distributions proposed by Biiyukkilig et al jl2| without any approximation. 
On the basis of the factorization theorem discussed above, the derivation of this kind of distribution is straightforward. 



III. COMPLETE DISTRIBUTION FORMALISM 



Let us begin with the complete distribution (Trp = 1) and the normalized average x = Trpx. This formalism 
of NSM is first proposed by Tsallis and has received little attention due to some problems (e.g. missing 
Legendre transformation and zeroth law, thermodynamic stability problem discussed by assuming extensive energy 
p3fl). Since several years, NSM has much evolved. Legendre transformation and the zeroth law in this formalism can 
be established by using the method of the references |l7J ■ The stability problem must be revisited with nonextensive 
energy satisfying Eq.(nh. So this standard formalism of probability theory remains a possible choices of NSM and 
still deserves to be studied. The reader will see below that the formalisms with normalized average can give the very 
distribution functions of Buyukkilic, et al. 

For nonextensive quantum gas, we have |?],[l^,^2| : 

p = i [1 - (q - 1)(3(H pN)} ^ = ie-^' (4) 
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where h! = ln [ 1 +( 1 ^g)^(^ ^ Af )l can j-, e called "deformed hamiltonian" . The grand partition function Z is then given 
by 



Z = lT[l-(q-l)/3{H-nN)]<=T (5) 
= Tr 

From the factorization theorem, we can write 



Z = Tre HS £» fc « (6) 

k rife 

where e& is the eigenvalue of h[ — lHl+i^|j^^ii£)] ; the deformed one-particle hamiltonian satisfying h' = h' i7 
and rife the occupation number of the one-particle state k. For boson and fermion, we obtain, respectively, 

oo 1 
fc rifc=0 fe k n k =0 k 

Then, it is straightforward to show that, just as in the conventional quantum statistics , 
_ l dQnZ) 1 1 

n ; = Trpni = - — — = = ; (8) 

(3 dei eP" ±1 [1 + (1 - q)/3(ei - fi)} ~ ± 1 

where e/ is the eigenvalue of the one-particle hamiltonian Hi. "+" and "-" correspond to fermion and boson, respec- 
tively. This result is just that given by Biiyukkilic et al JT3] with only a change (q — 1) — > (1 — q). 

IV. INCOMPLETE DISTRIBUTION FORMALISM 

For incomplete distribution with the ^-normalization, the Grand partition function Z is given by 

Z q = Tt[1- (l-q)/3(H - fxN)]^ (9) 
= Tre - * 3 ' 1 ' 

— e -qn k f3e k 
k n k 

— J^J e -q0n k e k 

k n k 

The deformed hamiltonians /i' and h[ are the same functions as given above with the transform (q — 1) — ► (1 — g) [ p"7| . 
The same machinery leads to 

ni — 1rp q ni — — 5 — = 5 . (10) 

e^±l [l + (g-l)/3( e/ - M )]i^±l 

These distributions are equivalent to Eq.(||) with only a difference in the g-dependence. Eq.(||) and Eq.([l0|) represent 
standard bosons and fcrmions satisfying Pauli exclusion principle. For example, ni < 1 for fcrmions. It should be 
noticed that, in this case, the Fermi surface at T = or (3 = oo is the same as in the conventional case and independent 
of the parameter q. 
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V. COMPLETE DISTRIBUTION WITH UNNORMALIZED AVERAGE 



On the other hand, in the formalism of complete distribution with unnormalized average, the things will be different 
because this average, unlike the normalized ones, can lead to nonadditive deformed internal energy even with additive 
deformed hamiltonian h ' . This formalism received a lot of attention and finally was found to show some puzzling 
properties 0. But up to now, there is no really solid reason for abandoning it definitely. We still see various kind 
of unnormalized expectations widely used in many fields, especially in financial problems |6[|24|| . It is also found that 
this formalism is the only one of NSM that can give standard statistical interpretation of heat and work j2j| . So it 
is still of interest to study this formalism in the circumstance of quantum systems. 

The partition function is still given by Eq.(||) with a transform (q — 1) — > (1 — q) [jf). The occupation number is 
calculated as follows 



ni = TT P q ni (11) 



Trn ; i-[l + (g-l)/3(ff-/xiV)]T^ 



k n k 

-1 () -i - . - D -qn k f3e k 



e -qn k /3t k \_ 



Zi ei^i - 1 

k n k 

Trp q 



g/3ei _ 1 



1 _T\(±~ e ~ 0ek ^ 
i ^ i 



e q0ei _ I 11 I _ e -9/3e 



for bosons and 



ni 



1 yr l+e^^ efc 

e qf3ei + 1 11 ( 1 + e -/3e fe )g C ) 



for fermions. Eq.(|Tl|) and (12) can be written as 

Q± Q± 

n l — a = 3 

e^ e i±l [l + ( 9 -l)/?( ei ±i 

where Q + = Trp 9 = J| fe ! ^ e e -'f,^ q ( or Q- = Ilfc ^Yzf^i^r") can be regarded as a parameter depending on q. Q± > 1 
Q± = 1 and Q± < 1 for q < 1, q = 1 and g > 1, respectively. These distributions seems interesting because they 
allow intermediate occupation number between that of bosons and fermions. In particular, at absolute zero, for 
"fermionlike" particle with "+" in Eq. (|l3|) , hi — Q when e; < fi and hi = when e/ > /i. This means that it is 
possible for several "fermions" to occupy an one-particle quantum state if Q > 1 or q < 1. Consequently, the Fermi 
surface ep at T = or (3 = oo changes as a function of Q or of the interaction between the particles : ep = 
where Ef is the conventional Fermi energy at q = 1 or Q = 1. This result can be compared to that of the fractional 
exclusion statistics Q for intermediate particles different from bosons and fermions. It is not surprising to see that 
nonextensive statistics has similar effect to that of fractional statistics describing interacting particle or elementary 
excitation |^,|l^], because the g-distribution is nothing but a result of long range interactions. However, the fact 
that only the NSM formalism with unnormalized average can give the intermediate quantum distributions seems to 
deserve further investigation. Quite interesting efforts have been made by some authors @,^6| to relate nonextensive 
statistics and the quantum distributions given by Eq.(||) or ([h]) to fractional exclusion statistics. 

It is obvious that if we assume the g-normalization for incomplete distribution, Q — 1 and Eq.(|T3|) will become 
Eq.©. 
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VI. CONCLUSION 



Summing up, on the basis of the factorization theorem of compound probability prescribed by thermodynamic 
equilibrium, the quantum gas distributions within different formalisms of NSM are investigated. It is shown that 
only the formalisms with normalized average can give distribution functions for standard bosons and fermions. The 
formalism with unnormalized average leads to a intermediate distribution similar to that of fractional exclusion 
statistics, with Fermi surface at T = depending on the parameter q. This different quantum properties of NSM due 
to normalizations remains something to be understood. A detailed study of the relation between the unnormalized 
quantum distribution and the fractional exclusion statistics would be of interest. 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 
[7] 

[8] 
[9] 
[10] 

[11] 
[12] 



[13] 
[14] 

[15] 



[16] 
[17] 

[18] 
[19] 
[20] 
[21] 

[22] 
[23] 
[24] 
[25] 

[26] 



Biiyiikkilig F., Demirhan D. and Giileg A., Phys. Lett. A, 
Physica A, 1997;240:657; Tirnakli U., Biiyiikkilig F. and 



Kroger H, Fractal geometry in quantum mechanics field theory and spin systems, Physics Report, 2000;323(2):82-181 
El Naschie MS, On the elctroweak mising angle in E°° Chaos, Solitons & Fractals, 2000;11(11):1803-1807 
Nottale L, Scale relativity and fractal space-time Chaos, Solitons & Fractals, 1966;7(6):877-938 

Ahmed E, Hegazi AS, On infinitesimally deformed algebra and fractal space-time theory Chaos, Solitons & Fractals, 
2000; 11 (11): 1759 

Mandelbrot B., Fractals in Finance, Springer, Berlin, 1997 

C. Tsallis, Chaos, Solitons & Fractals, 1995;6:539. 

Tsallis C, J. Statis. Phys., 1988;52:479; Curado E.M.F., and Tsallis O, J. Phys. A: Math. Gen. 1991;24:L69; Tsallis C, 
Mendes R.S. and Plastino A.R., Physica A, 1998;261:534 
Badiali J.P., Phys. Rev. E, 1999;60:2533 
Wu Y.S., Phys. Rev. Lett, 1994;73:922 
Haldane F.D.M., Phys. Rev. Lett, 1991;67:937 
Monteiro MR, Roditi I., Mod. Phys. Lett. B, 1993;7:1843 
Biiyiikkilig F. and Demirhan D., Phys. Lett. A, 1993;181:24 
1995;197:209; Tirnakli U., Biiyiikkilig F. and Demirhan D. 
Demirhan D., Phys. Lett. A , 1998;245:62; Torres D.F. and Tirnakli U., Physica A, 1998;261:499; Tirnakli U. and Torres 

D. F., Physica A, 1999;268:152 

Biiyiikkilig F. and Demirhan D., Euro. Phys. J. B, 2000;14:705 

Kaniadakis G., Lavagno A. and Quarati P., Modern Phys. Lett. B, 1996;10:497; Kaniadakis G., Lavagno A. and Quarati 
P., Nucl. Phys. A, 2001;680:94c-97c 

Tsalli s C., Entropic non extensivity : a possible measure of complexity, Chaos, Solitons & Fractals, 2001;13:371 and e- 
print, cond- mat/0010150 Tsallis C, Nonextensive statistics : theoretical, experimental and computational evidences and 
connections, Brazilian J. Physics, 1999;29:1 
Wang Q.A. and Le Mehaute A., Phys. Lett. A., 1998;242:301 

Wang Q.A., Pezeril M., Nivanen L., Le Mehaute A., Chaos, So litons fc Fractals, 20 02;13:131, e-print : :ond-mat/0010294 ; 

Wang Q.A., Chaos, Solitons & Fractals, 2001;12:1431, e-print, |cond-mat/0009343 

Lenzi E.K., Mendes R.S., da Silva L.R. and Malacarne L.C., Physics Lett. A., 2001;289:44 

dos Santos R.J.V, J. Math. Phys., 1997;38:4104; Abe S., Physics Lett. A., 2000;271:74 

Nakamichi A., Joichi I., Iguchi O. and Morikava M., Chaos, Solitons & Fractals, 2002;13:595 

S. Abe, Phys. 7?et;.£:,63(2001)061105; Qiuping A. Wang, Laurent Nivanen, Alain Le Mehaute and Michel Pezeril, Note on 
Abe's general pseudoadditivity for nonextensive systems, e-print : cond-mat/01 11541 (submitted) 
Curilef S., Z. Phys. B, 1996;100:433 

Ramshaw John D., Phys. Lett. A., 1995;198:119; Guer beroff G.R. and Ra ggio G.A., Phys. Lett. A., 1998;214:313 



C. Anteneodo, C. Tsallis and A.S. Martinez, e-print : :ond- mat/0109203 

Wang Q.A., Unnormalized nonextensive expectation value and zeroth law of thermodynamics, e-print : cond-mat/0111238 
(submitted) 

Biiyiikkilig F., H. Uncu and D. Demirhan, Thermodynamics of the g-on gas, eprint: cond-mat/0103560 



5 



